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ABSTRACT 

Most black hole binaries show large changes in X-ray luminosity caused primarily by variations in mass ac- 
cretion rate. An important question for understanding black hole accretion and jet production is whether the 
inner edge of the accretion disk recedes at low accretion rate. Measurements of the location of the inner edge 
(R m ) can be made using iron emission lines that arise due to fluorescence of iron in the disk, and these indicate 
that Ri n is very close to the black hole at high and moderate luminosities (> 1% of the Eddington luminosity, 
^Edd)- Here, we report on X-ray observations of the black hole GX 339—4 in the hard state by Suzaku and the 
Rossi X-ray Timing Explorer (RXTE) that extend iron line studies to 0.14% Ledd and show that R[ n increases 
by a factor of >27 over the value found when GX 339-4 was bright. The exact value of /?;„ depends on the 
inclination of the inner disk (i), and we derive 90% confidence limits of /?;„ > 35/?,, at i = 0° and R m > \15R g 
at i = 30°. This provides direct evidence that the inner portion of the disk is not present at low luminosity, 
allowing for the possibility that the inner disk is replaced by advection- or magnetically-dominated accretion 
flows. 

Subject headings: accretion, accretion disks — black hole physics — stars: individual (GX 339^4) — X-rays: 
stars — X-rays: general 



1. INTRODUCTION 

Accreting black holes exhibi t a variety of spectral states 
dMcClintock & R emillard 2006), most often showing strong 
thermal emission from a n optically thick accretion disk 
(Shakura & Sunyaevlll973l) when they are at their brightest. 
As luminosities drop below a few percent of the Eddington 
value (LEddX the level of thermal emission decreases, and 
the spectrum harden s, signaling a transition to the hard state 
dKalemci et al.l l2004). This state is also characterized by hav- 
ing a high level of X-ray variability and steady, self-absorbed 
compact jets, which are often observed in the radio band 
dFenderll2001h . 

In theory, it is predicted that the cool gas from the 
accreti on disk will evaporate as the mass accretion rate 
drops (Me yer. Liu & Meyer-Hofmeisterll20"00l) . The implied 
increase in the inner radius of the optically thick disk, 
R ln , is critical to the idea that a quasi-spherical advection- 
domin ated accretion flow ( ADAF) forms around the black 
hole dNaravan & Yil |1994|) . This model requires that R in 
increases to ~ 100- 1000 R g (lEsin. McClintock & Naravanl 
1997), where the gravitational radius is R g = GM/c 2 , G and 
c are constants, and M is the black hole mass. 

Recently, observational efforts to test the truncated disk 
picture have used two techniques: Modeling the thermal 
emission from the disk and measuring the iron Kg emis- 
sion line in the reflection component (Light man & Whitel 
1988). The for mer technique is challenging because of un- 
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certainties arising from electron scattering, irradiation of the 
disk, and the torque boundary condition at the disk's in- 
ner edge as well as the difficulty of modeling a soft X- 
ray component that is strongly affected by interstellar ab- 
sorption. Modeling the disk's thermal emission to estimate 
the inner radius (Rm) has le d to mixed results with some in 
favor of an increas e in R{ n oGierlinski. Done & Page! 120081: 
I Cabanac et al.l200"9[) and others in d isagreement dRvkoff et al.l 
l2007t iReis. Miller & Fabianl 12009). The iron line provides a 
measurement of R[ n because it is Doppler-broadened due to 
relativistic motion of the accretion disk material and gravita- 
tionally redshifted due to the black hole's gravitational field. 
Broad iron lines have been seen from both ste l lar mass and 
super massive black holes (Tanak aet al.l 119951: iFabian et al.l 
2009). One challenge to the relativistic interpretation in the 
case of stellar mass black holes is tha t the iron line shape 
may be affected by scattering in a wind dLaurent & Titarchukl 
120071: iTitarchuk. Laurent & Shaposhnikovll2009i) . However, 
argumen ts in favor of the relativistic interpretation are pre- 
sented in iMillerl d2007l) . 

For stellar mass black holes at high luminosities, the iron 
line profiles have been used extensively, and values of Ri n 
less than the radius of the innermost stable circular orbit for 
a non-rotating black hole, Rj n < 6 R„, have been measured 
for several systems dMiller et alj2002tl 2"b04. 2008), including 
GX 3 3 9—4 . The broad iron line has, somewhat surprisingly, 
persis ted into the hard state dMiller et al .1120061: iTomsick et al.l 
l2008l:lHiemstra et al.ll2009l) . However, for GX 339-4, the pre- 
vious detections of the iron line in the hard state were obtained 
at luminosities above 1% LEdd, and it remains to be seen if 
there is evolution in the iron line profile at lower luminosities. 

GX 339^4 is one of the most active black hole tran- 
sients known, with 4 major outbursts in the past 7 yr . 
The distance to GX 339-4 is -8 kpc dHynes et all 12004 
IZdziarski et al.l 12004]). and the orbital peri od is 1.7 days 
dHynes et all l2003b Levine & Corbetl 120061) The lower 
limit on the black hole mass is >6M^ dHvnesetal.l l200l 
iMufioz-Darias. Casares & Martfnez-Paisl 2008), and here, we 
adopt a value of 10M©, corresponding to Lndd ~ L3 x 10 39 
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ergs s -1 . GX 339—4 has allowed for some of the most in-dept h 
studies of compact jets (ICorbel et al.l2000tlCoriat et alJ2009h . 
Here, we use Suzaku and RXTE observations of GX 339-4 to 
study the iron line profile at 0.14% LecM- 

2. OBSERVATIONS 

We observed GX 339^4 in the X-ray band over a period 
of 2.2 days during 2 008 Septemb er 24-27 with the Suzaku 
(Mits uda etaTll2007l) and RXTE dBradt, Rothschild & Swankl 
11991 " satellites. The observation was made 1.6 yr after the 
peak of its 2007 outburst during a time when the source was 
active but faint dRussell et al.l 2008; Kong 2008). We ana- 
lyzed the data using the HEASOFT v6.7 software, and the 
most recent instrument calibrations as of 2009 September 
25 for Suzaku and 2009 August 19 for RXTE. Suzaku ob- 
served GX 339^4 continuously over the 2.2 day period (ex- 
cept for Earth occulations), and we accumulated 105,000 s 
of on-source time for the X-ray Imaging Spectrometer (XIS) 
and 106,800 s for the PIN layer of the Hard X-ray Detector 
(HXD). The Suzaku data are all contained in observation ID 
403067010. RXTE made five shorter (2,300 to 3,500 s) obser- 
vations during the same time period. The observation IDs are 
93702-04-01-01, 93702-04-01-02, 93702-04-01-03, 93702- 
04-02-00, and 93702-04-02-01, and we obtained a total Pro- 
portional Counter Array (PCA) exposure time of 14,640 s and 
a total High-Energy X-ray Timing Experiment (HEXTE) ex- 
posure time of 4,625 s. We also observed GX 339^4 in the 
radio band at the Australia Telescope Compact Array (ATCA) 
on 2008 August 18 and September 28. 

For observations of relatively faint sources close to the 
Galactic plane, an additional step is required for the PCA 
spectrum to subtra ct any other emission in its field-of-view 
( Jahoda et al. 2006). To do this, we obtained data from RXTE 
observations of GX 339-4 whe n the source was in quiescence 
dGallo. Fender & Corbel 2003). These observations occurred 
between 2001 March and 2002 February and also on 2003 
September 29. The average quiescent count rate is ~6% of 
that measured during the 2008 September observations. We 
subtracted the quiescent PCA spectrum for the spectral fits 
below. 

3. RESULTS 

3.1. General Properties and Spectral State Identification 

We first used the X-ray and radio properties to iden- 
tify the spectral state. We produced the energy spec- 
trum shown in Figure QJ, and fitted it with an absorbed 
power-law model. To model the absorption, we used 
elemental abundances that ap proximate interstellar values 
dWilms. Allen & McCravl l2000h and atomic c ross sections 
from iBalucinska-Church & McCa mrnonl (1 19921) . For the XIS 
spectra, we did not include the data in the 1.7-1.9 keV range 
due to uncertainties in the instrument calibration near the Sili- 
con K-edge. The overall normalizations between instruments 
were left as free parameters in the model, and the fluxes given 
below are based on those measured by XIS0+XIS3. The mea- 
sured model parameters are the interstellar column density, 
Af H = (6.79 ± 0.07) x 10 21 cm" 2 , the power-law photon index, 
T = 1.573 ±0.006, and the unabsorbed 0.4-12 keV flux of 
(9.08 ±0.03) x 10" 11 ergs cm" 2 s" 1 (90% confidence uncer- 
tainties, A\ 2 = 2.7). Although this model provides a good 
description of the continuum, the quality of the fit is not good 
(X 2 /v = 1044/597) due to the presence of a strong iron Ka 
line near 6.4 keV. 




1 10 100 

Energy (keV) 



FIG. 1 . — The Suzaku and RXTE energy spectrum for the entire duration of 
the 2008 September 24-27 observation of GX 339—4. (a) The counts spectrum 
fitted with an absorbed power-law model, (b) The residuals for the power-law 
fit to the spectrum given in terms of the contribution from each data point to 
the goodness-of-fit parameter, \. The largest residuals are seen near 6-7 keV 
and are due to the iron Ka emission line. 



To characterize the X-ray timing properties, we produced 
XIS light curves in the 0.4-2.5 keV and 2.5-12 keV energy 
bands with 8 s time resolution and PCA light curves in the 3- 
25 keV band with 16 s resolution. An inspection of the light 
curves shows that while they show strong variability and flar- 
ing, they are very similar in the different energy bands, indi- 
cating that the spectrum changes very little with flux. We also 
produced a PCA power spectrum at frequencies from 0.005 to 
64 Hz. The power spectrum is well-described by a power-law 
with an index of 1.11 ±0.03, and a fractional rms value of 
55% ±2% (0.01-10 Hz). 

The ATCA radio observations show flux densities of 1 . 10 ± 
0.10 and 1.18 ±0.10 mJy at frequencies (v) of 4.80 and 
8.64 GHz, respectively, on August 18. For the same fre- 
quencies, the September 28 flux densities are 0.9 ± 0.2 and 
1.1 ±0.2 mJy. These measurements indicate spectral indices 
of a = 0.13 ±0.22 and a = 0.3 ±0.5, where the flux density 
is given by S„ oc v a . Self-absorption of the compact jets that 
are seen in the hard state causes flat or slightly inverted radio 
spectra (a > 0.0), and our measurements are consistent with 
this. 

Thus, the hard energy spectrum, the high level of timing 
noise, and the evidence for the presence of a compact jet are 
all consistent with the source being in the hard state during 
our observation. 

3.2. The Iron Emission Line and the Inner Radius of the Disk 

The residuals shown in Figure [Tp indicate the deviation of 
the data from a power-law, and the presence of an iron Ka 
emission line is clear. When we add a Gaussian to model the 
emission line, we measure a line energy of 6.45^' keV, con- 
sistent with neutral to moderately ionized iron, a line width of 
a = 0.14^$* keV, and an equivalent width (EW) of 77+}g eV. 
With the addition of the iron line, the quality of the fit im- 
proves dramatically to \ 2 l v = 792/594. The total 1-100 keV 
unabsorbed flux is 2.4 x 10~ 10 ergs cm" 2 s" 1 , which is a fac- 
tor ~9 lower than the lowest level at which an iron line was 
previously detected for GX 339-4 dTomsick et alj|2008l) . and 
this flux corresponds to a luminosity of 0. 14% Ledd- 
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Before discussing the implications of the presence of this 
narrow iron line, it is critical to determine if the iron line 
could be related either to poor background subtraction or to 
emission from other sources in the Galactic plane. For the 
XIS detectors, we examined the background spectrum from 
two 4'. 4 x 3'. 7 rectangular regions on the detectors and veri- 
fied that only internal ba ckground lines of the detector appear 
(see lKovama et al. 2007, for information about the XIS inter- 
nal background) without any evidence for other background 
lines with strong emission in the iron Ka region. 

Emission from the Galactic ridge includes iron Ka emis- 
sion with the m ost prominent line being due to He-like 
iron at 6.7 keV dKovama et al l Il986t iKaneda etalJ 11997b 
iRevnivtsev et al.ll2009t) . Although the emission is very strong 
in the Galactic center region, it decreases for lines-of-sight 
away from the Galactic center and drops especially rapidly 
with Galactic latitude (b). In the Scutum region, at a Galac- 
tic longitude of / = 2 8.5°, a scale height of 0.5° is estimated 
(IKaneda et al.1 119971) . Thus, at the position of GX 339^1 
(I = 338.9°, b = -4.3°), the Galactic ridge emission is expected 
to be weak, but possibly not negligible. 

The fact that the XIS0+XIS3 background spectrum does not 
show an emission line at 6.7 keV suggests that we are not de- 
tecting Galactic ridge emission in our observation of GX 339- 
4. As an additional check, we produced a new GX 339-4 
spectrum using a circular extraction region centered on the 
source with a radius of 0'.86, which is five times smaller than 
the 4'. 3 radius region used previously. This causes a reduc- 
tion in the source count rate by a factor of 2.1 while reduc- 
ing the background by a factor of 25. Thus, the strength of 
any background features will decrease by an order of mag- 
nitude. Fitting the new GX 339^1 XIS spectrum with an 
absorbed power-law model and inspecting the residuals still 
clearly shows an iron line at 6.4 keV. Adding a Gaussian to fit 
the iron line, we measure an EW of 73t|J eV, which is consis- 
tent with no change from the EW of 71+jo eV that we measure 
with XIS with the larger extraction region. Based on this and 
the background spectrum discussed above, we conclude that 
the iron line is from GX 339-4. 

To use the shape of the iron line to constrain R[ n , we return 
to fitting the full spectrum shown in Figure[T| and we replace d 
the Gaussian component with the laor model (lLaorll 1991b . 
which accounts for the relativistic effects near a rotating black 
hole. The line energy and EW are 6.47^Q3 keV and 72+' eV, 
respectively. The other model parameters are /?; n , the incli- 
nation of the inner disk, i, and the emissivity index, q, which 
is a power-law index that sets how the line-emissivity (J) of 
the disk changes with radius according to J oc r~ q . Although a 
broad line allows for good constraints on all 3 parameters, this 
is not the case for a narrow line. Thus, we restricted the range 
of q to be between 2, which corresponds to a value where 
the co ntribution to the iron line from large radii begins to di- 
verge (Laor 1991), and 3, which is consistent with the value 
obtained from previo us measurements of the GX 339-4 iron 
line in the hard state dMiller et alj|2006l 120081: iTomsick et all 
120081: [Reis et al. 20081). The inclina tion has been prev iously 
measured to be i = 18 ±2 degrees dMiller et al.ll2008l) . For 
Ri n , the spectral fits indicate 68% and 90% confidence lower 
limits of >84 R g and >65 R g , respectively, if z" = 18°. As a 
value of Rj n = 2.4 R e was obtained when the source was bright 
dMiller et a l. 2008), our results indicate that the inner radius 
changes by a factor of >27. Figure |2]illustrates the huge dif- 
ference between the profile that we measure and a profile with 
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FIG. 2. — Profile of the iron Ka emission line from GX 339^4 measured 
by XIS0+XIS3 (in grey). The profile is shown in terms of a data-to-model 
ratio, where the model is an absorbed power-law fitted to the Suzaku and 
RXTE spectra, excluding the 4.0-8.0 keV region. We determined the model 
profiles (in black) by restoring the 4.0-8.0 keV data and fitting the entire 
spectrum with a model consisting of an absorbed power-law plus a relativistic 
laor iron line. For the profiles shown, we used a disk inclination of 18° 
and an emissivity index of q = 3. The solid black line shows a profile with 
Rin = 105 R g , which corresponds to the best fit value. The other two profiles 
(dashed and dotted black lines) were obtained by fixing Rj n to 30 R g and 
2.4 R g , respectively, and then refitting the spectrum, allowing the line energy 
and normalization to adjust. 



R m = 2AR g . 

Fixing the disk inclination to i= 18° is appropriate for con- 
straining the ratio of inner radii at 0.14% and 12-15% Ledd 
since we do not expect the disk inclination to change signifi- 
cantly with luminosity, being set either by the binary inclina- 
tion or the spin axis of the black hole. However, for obtaining 
a physical value of Ri n , we need to consider the inclination. 
We refitted the full spectrum with the power-law plus laor 
model, allowing i to be a free parameter and q to be free in 
the range 2-3. Then, we performed a grid search of incli- 
nations covering 0° to 40° and values of R m covering 10 to 
210 R g , Figure [3] shows the results in terms of the 68% and 
90% confidence contours. The contours show that at lower 
inclinations, the constraint on R m becomes somewhat weaker, 
but even at i = 0°, the results indicate a truncated disk with 
R\ n > 35R g (90% confidence). On the other hand, at disk in- 
clinations above 18°, the inferred values of R[ n rise rapidly. 
For example, at i = 30°, the limit is R m > H5R g . 

There are at least two reasons why our limits on R m are 
conservative in the sense that the di sk may be even more trun- 
cated. First, Ross & Fabian (2005) show that, especially for 
the disks around stellar mass black holes, thermal motions 
of the disk material can lead to significant broadening of the 
iron line. If this is the case for GX 339^1, then less Doppler 
broadening from bulk motion in the accretion disk would be 
required. Second, it has been pointed out that if the disk in- 
clination is really as low as 18° and the disk inclination is the 
same as the binary inc lination, then this wo uld imply a black 
hole mass of 200 M (Cabanac et al. 2009), which is proba- 
bly unrealistically high. Thus, either the disk is warped or the 
inclination is significantly higher than 18°. 
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FIG. 3. — Confidence regions for two laor iron line parameters: The disk 
inclination and the inner radius of the disk. The dark grey region shows the 
68% confidence (Ax 2 = 2.30) error region and the light grey region shows 
the 90% confidence (Ax 2 = 4.61) error region. 



4. DISCUSSION 

These results provide the most direct and quantitative evi- 
dence to date for the truncation of the accretion disk for stellar 
mass black holes in the hard state at low luminosities. Fig- 
ure^ compares our constraint on R in at 0. 14% Lndd to previ- 
ous measurements at higher luminosities. The data show that 
the inner edge of the disk moves sharply outward as the lumi- 
nosity decreases from 1% to 0.1% I^dd- In addition, the drop 
in the iron line EW shown in Figure |4j? is consistent with the 
increase in since one expects the line to become weaker 
if the disk is truncated. The EW evolution and the fact that 
the line is well-described by a single component provide sup- 
port for the interpretation that the narrow line comes from the 
disk. Considering the broad lines seen in the brighter part of 
the hard state and the disk truncation that we see here, the 
overall evolution is very similar to the theoretical prediction 
that the disk evaporation will start in the middle of the disk 
and that the entire inner di sk will evaporate when the lumi- 
nosity reaches ~0.1% Z^dd dTaam et al.ll2.Q08l) . 

While the study of stellar mass black holes, such as 
GX 339-4, provides a detailed look at the evolution of R in , 
there is iron line evidence for truncated disks around super- 
massive bl ack holes in Active G alactic Nuclei (AGN). For 
NGC 4258. Reynolds et al.l(l2009l) measure a narrow line with 
an EW of 45 eV at a luminosity of 10~ 3 % Z^dd and infer that 
Ri n > 3 x 10 3 R g . There is also evidence for a truncated disk in 
NGC 4593, and for this AGN, there is possible evidence for 
a cha nge in Ri a over a period of 5 yr (Markowi tz & Reevesl 
120091) . 

Our finding for GX 339^- is also important because a large 
increase in R[ n is required for the presence of an ADAF, and 
our result makes the ADAF model viable for the fainter part 
of the hard state. However, it is known that the ADAF model 
does not give a complete physical description of the system 
since it does not incorporate the compact jet. These jets have 
now been detected from GX 339^- in the hard state both 
when the disk is truncated (this work) and when the disk is 



not truncated dMiller et al.l l2006t iTbmsick et al.l l2008). The 
possibility that disk truncation occurs at a different mass ac- 
cretion rate than the transition to the hard state and the turn- 
on of the compact jet has not been considered in most the- 
oretical work. This finding strongly constrains models such 
as the magn etically-dominated accretion flow (MDAF) model 
dMeierll2005l) where jet production and properties depend on 
the accretion geometry. 
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FIG. 4. — Iron line measurements for GX 33<M- over a 0.5-100 keV 
luminosity range from 0.14% to 15% LEdd- (a) The constraints on the in- 
ner radius of the accretion disk, including all GX 339^1 observations made 
with instruments with energy resolution better than 200 eV for which a de- 
termination of i?j„ was obtained using relativistic reflection models. The 
four measureme nts of R m in the ra nge 2-5 R ? bet ween 1.3% Ledd and 15% 
L Edd use Suzaku Miller et all2008r>. XMM-Newt on Miller et alJ200ll2006l: 
IReis et 31112003) and Swift jTomsick et alj|2008l) . The grey 90% confidence 
upper limit is a Swift measurement where broad features in the reflection com- 
ponen t are detected, but the iron line is not clearly detected I Tom sick et all 
2008). The dashed line marks the central value obtained at 12-15% LEdd- 
The measurements indicate that a large increase in R m occurs between 0.5% 
LEdd and 0.14% LEdd- (b) The equivalent width of the iron line vs. luminosity. 
The value we obtain at 0.14% LEdd is the lowest value, which is consistent 
with a drop in the strength of the reflection component due to truncation of 
the disk. 
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